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Abstract: Chronic obstructive pulmonary disease (COPD) is a major global health problem. 
It results from chronic inflammation and causes irreversible airway damage. Levels of different 
serum cytokines could be surrogate biomarkers for inflammation and lung function in COPD. 
We aimed to determine the serum levels of different biomarkers in COPD patients, the asso- 
ciation between cytokine levels and various prognostic parameters, and the key pathways/ 
networks involved in stable COPD. In this study, serum levels of 48 cytokines were examined 
by multiplex assays in 30 subjects (control, n=9; COPD, n=21). Relationships between serum 
biomarkers and forced expiratory volume in 1 second, peak oxygen uptake, body mass index, 
dyspnea score, and smoking were assessed. Enrichment pathways and network analyses were 
implemented, using a list of cytokines showing differential expression between healthy controls 
and patients with COPD by Cytoscape and GeneGo Metacore™ software (Thomson-Reuters 
Corporation, New York, NY, USA). Concentrations of cutaneous T-cell attracting chemokine, 
eotaxin, hepatocyte growth factor, interleukin 6 (IL-6), IL-16, and stem cell factor are signifi- 
cantly higher in COPD patients compared with in control patients. Notably, this study identifies 
stem cell factor as a biomarker for COPD. Multiple regression analysis predicts that cutaneous 
T-cell-attracting chemokine, eotaxin, IL-6, and stem cell factor are inversely associated with 
forced expiratory volume in 1 second and peak oxygen uptake change, whereas smoking is 
related to eotaxin and hepatocyte growth factor changes. Enrichment pathways and network 
analyses reveal the potential involvement of specific inflammatory and immune process pathways 
in COPD. Identified network interaction and regulation of different cytokines would pave the 
way for deeper insight into mechanisms of the disease process. 
Keywords: COPD, Bio-Plex assay, biomarkers, pathways, networking 

Introduction 

Chronic obstructive pulmonary disease (COPD) is a major global health problem. 
In 2002, it was the fourth worldwide leading cause of death, and it is anticipated 
to be the third leading cause of death by 2030. 1 COPD results from chronic inflam- 
mation and eventual irreversible damage to the airways. Several factors, including 
exposure to environmental pollution and smoking, contribute to the airway damage 
and obstruction. 2 Inflammatory response in COPD is characterized by the increased 
number of macrophages, neutrophils, and cytotoxic T lymphocytes in the airways and 
lung parenchyma. Airway inflammation further results in systemic inflammation and 
other COPD-related manifestations. 3 3 
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Spirometry (forced expiratory volume in 1 second; 
FEVj), exercise tests (peak oxygen uptake; V0 2 ), and dys- 
pnea scores (a scale to determine the degree of difficulty 
in breathing) are used for diagnosing and determining 
the prognostic factors of COPD. These tests are also used 
for evaluating the treatment efficacy in COPD patients. 
However, the test results do not fully reflect the burden of 
COPD. 6 7 Aaron et al have studied 19 serum markers and 
found that only C-reactive protein, myeloperoxidase, and 
vascular endothelial growth factor show intra- and inter- 
patient reliability and correlation with disease severity. 8 
Franciosi et al performed a meta-analysis by considering 
the spirometric, demographic, clinical, cytological, and 
biochemical variables in stable COPD patients. Their meta- 
analysis indicates that the concentrations of tumor necrosis 
factor (TNF) and C-reactive protein have suggestive cor- 
relation with disease severity. Further, the existence of high 
variability among biomarker levels in different studies has 
also been observed. 9 Therefore, additional studies with a 
well-defined group of patients are needed to identify the sets 
of potential biomarkers and pathways involved in COPD 
from different populations. 

Cytokines play a critical role in COPD-associated 
inflammation. 10 Interactions among these cytokines are 
complex, as some cytokines can have different actions under 
a specific context. 11 Thus, studying a panel of cytokine 
biomarkers, their correlation with pulmonary physiologi- 
cal functions and systemic parameters and the involvement 
of various pathways and networks would provide a better 
understanding of the mechanisms underlying COPD. This 
will also help to identify a nonredundant set of biomarkers for 
COPD. Measuring cytokine panels from the same samples in 
multiplex assays will be advantageous for reducing technical 
variability. Only a limited number of studies have used panels 
of multiple serum cytokines to determine biomarkers in 
COPD patients simultaneously. 812 Furthermore, no studies 
with multiple biomarkers and pathway analysis have been 
reported on COPD patients of Indian descent. 

In the present study, we measured the serum levels 
of 48 cytokines in stable COPD and healthy subjects by 
immunobead Bio-Plex assays (Bio-Rad Laboratories Inc., 
Hercules, CA, USA) and analyzed their association with 
FEVj (percentage predicted), maximum oxygen uptake 
(peak V0 2 predicted), smoking index, dyspnea scores, and 
body mass index (BMI) by multiple regression analysis. 
Furthermore, our study identifies that changes in a few key 
cytokine concentrations explain the changes in lung func- 
tion related parameters (FEVj and peak V0 2 and smoking 



index). We also assessed some of the key biomarkers and 
their common regulatory factors involved in inflammatory 
and immunological pathways in COPD. 

Materials and methods 

Ethics statement 

This study was approved by the Ethics Committee of All 
India Institute of Medical Sciences, New Delhi, India. Prior 
informed written consent was obtained from all participating 
healthy and COPD volunteers. Blood work and spirometric 
tests were conducted under the guidelines of the ethics 
committee. 

Patient selection 

COPD patients were diagnosed according to the guidelines 
of Global Initiative for Chronic Obstructive Lung Disease 
(GOLD) criteria. 13 They were either current smokers with 
COPD or ex-smokers with COPD. Individuals older than 
40 years were recruited for this study and were categorized 
into two main groups: Group 1 included healthy controls 
with no evidence of COPD (n=9); four and five participants 
were smokers and nonsmokers, respectively, and group 2 
included patients diagnosed with COPD (n=21); eleven of 
whom were ex-smokers with COPD and ten of them were 
current smokers with COPD. All the COPD patients selected 
for the study were in stable condition (no exacerbation for at 
least 1 month). Participants were excluded specifically when 
they had a history of COPD exacerbation within the last 
1 -month period, asthma, or were taking oral steroids or had 
any conditions that were associated with inflammation such 
as infection, cancer, congestive heart failure, end-stage renal 
disease, or connective tissue disorders that could affect the 
levels of COPD biomarkers. 

Measurement of pulmonary functions 

FEV ( (percentage predicted) and FVC (forced vital capacity) 
were measured with a calibrated spirometer (Spiro Air, Medi- 
soft; PK Morgan Ltd., Kent, United Kingdom), using stan- 
dard methodology at the time of enrollment of the subjects. 14 
This pulmonary function test value is critical for diagnosing 
obstructive and restrictive lung diseases. 

Calculation of BMI 

BMI is an independent prognostic factor for disease severity 
and survival after the diagnosis of COPD. 15 BMI is a measure 
of human body shape/obesity based on an individual's weight 
and height. Each body mass index score was calculated by the 
following formula: BMI = weight (kg)/(height [m]) 2 . 
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Symptom-limited incremental 
exercise test 

In patients with COPD, a symptom-limited incremental step 
test elicited the maximum cardiopulmonary and metabolic 
responses. 16 A symptom-limited incremental exercise test was 
performed on a bicycle ergometer, using 10 W ramp protocol. 
Exercise was performed while seated on an electronically 
braked cycle ergometer (Corival; Lode, Groningen, the 
Netherlands). Heart rate was monitored continuously 
throughout the protocol using a heart rate monitor (POLAR 
[Polar Electro Oy, People's Republic of China]). The polar 
device used in this study had a T3 1 transmitter and a chest belt 
(N2965) to transmit heart rate to the system. Breath-by-breath 
measurements of oxygen consumption (V0 2 ) were made by 
custom software (vacuumed/vista MX [VacuMed, Ventura, 
CA, USA]). Standard 12-lead electrocardiograms were 
obtained at rest, during exercise and in recovery phase. Blood 
pressure was monitored using a standard cuff sphygmoma- 
nometer at rest, every 3 minutes during exercise, and every 
2 minutes during recovery. At the end of exercise, the reasons 
for termination of exercise were obtained from the patient. 
The VO, max was the highest V0 2 observed during exercise. 
The exercise was also terminated at exhaustion (intolerable 
dyspnea, as indicated by the patient); severe desaturation 
(peripheral capillary oxygen saturation <80%); demand by 
the patient because of leg cramps, chest pain, or discomfort 
or request by patients for any other reason. 

Evaluation of dyspnea 

Dyspnea is one of the common symptoms in COPD either 
at rest or under conditions of exercise, and this applies to 
all severities. Dyspnea severity was assessed and scored 
using the standard Modified Medical Research Council 
questionnaire. 1718 

Measurements of blood markers 

A volume of 5 mL of whole blood was collected into 
a plain tube and allowed to clot for 1 hour. The sera 
samples were aliquoted after centrifugation at 1,000 x g 
for 10 minutes and stored at -80°C until further analysis. 
Concentrations of serum cytokines were analyzed using 
two complementary Bio-Plex suspension array systems 
(Bio-Plex Pro Human Cytokine Group 27-Plex Panel and 
Bio-Plex Pro Human Cytokine Group 21-Plex Panel) to 
cover the range of all the cytokine biomarkers potentially 
involved in the pathophysiology of COPD. Forty-eight 
biomarkers were assessed simultaneously, using the Bio- 
Plex system. The selection of specific cytokines in the 



study was based on the previously available reports and 
the involvement of these cytokine biomarkers in COPD. 19 
Assays were performed in duplicate by following the stan- 
dard operating protocol provided by the Bio-Plex Multiplex 
cytokine assay. 

Briefly, anticytokine antibody-conjugated beads were 
added to individual wells of a 96-well filter plate and adhered 
using vacuum filtration. After washing, 50 |iL prediluted stan- 
dards and serum samples were added into respective wells, 
and the filter plates were shaken at 300 rpm for 30 minutes 
at room temperature. Thereafter, the filter plates were washed 
and 25 \iL prediluted multiplex biotin-conjugated detection 
antibody was added and incubated for another 30 minutes. 
After washing, 50 |iL prediluted streptavidin-conjugated PE 
was added for 1 0 minutes, followed by an additional wash and 
the addition of 125 jiL Bio-Plex assay buffer to each well. 
Then, filter plates were analyzed using the Bio-Plex Protein 
Array System, and concentrations of each cytokine were 
determined using software (Bio-Plex Manager version 6.0). 
Standard curves were generated for each biomarker. The line 
of best fit was determined for standard curves by standard 
recovery methods and by calculating the concentration of 
each standard. 20 

Statistical analysis 

Cytokines' levels in sera showed no differences between 
smokers and nonsmokers within the control group. There were 
also no statistical differences detected between ex-smokers 
with COPD and current smokers with COPD. Therefore, 
we compared the pooled samples (healthy controls versus 
COPD), using the two-tailed Mann- Whitney test, as outlined 
in the experimental design (Figure SI). The association 
of identified serum biomarkers with FEVj (percentage 
predicted), peak V0 2 (predicted), smoking index, dyspnea 
scores, and BMI were studied in all subjects, using Spear- 
man correlation analysis. In addition, a step-wise multiple 
regression analysis with weighted regression equations has 
also been used for calculating coefficients for FEV [ (percent- 
age predicted), peak V0 2 (predicted), smoking index, dyspnea 
scores, and BMI with serum biomarkers. Multiple regression 
analysis determines the relationship between a dependent or 
criterion variable of interest (Y) and a set of K independent 
variables or potential predictor variables (X x , X v X v ... , X^, 
where the scores on all variables are measured for n number 
of cases. A multiple regression equation for predicting Y is 
expressed as Y = A + (B^XJ + {B^C 2 ) + ... B k X k . A stepwise 
multiple regression analysis eliminated the parameters 
that provided insignificant contribution to the model. The 
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coefficient for each factor was obtained from the final model 
with the lowest P-value. Calculated and observed values 
were plotted in a regression line after calculating a common 
factor for each equation. All the statistical analyses were 
performed using the Prism version 5 (GraphPad Software, 
Inc., La Jolla, CA, USA) and SPSS (version 16; SPSS Inc., 
Chicago, IL, USA) software. 

Enriched pathways and network analysis 

The mean fold changes of the detected biomarkers were 
calculated between COPD and control subjects. Cytokines 
showing 1.5-fold changes were used for biological pro- 
cesses, pathways, and network enrichment analyses. The 
purpose of this analysis was to identify the enriched bio- 
logical processes and the pathways associated with COPD. 
Interaction networks among the significantly elevated 
cytokines in COPD were constructed using Cytoscape 
and GeneGo Metacore™ software (Thomson Reuters, 
St Joseph, MI, USA). 



Results 

Biomarkers CTACK, eotaxin, HGF, 
IL-6, IL- 1 6, and SCF are elevated 
in the sera of stable COPD patients 

We assessed 48 candidate biomarkers in the serum samples 
of healthy control and COPD subjects by Bio-Plex assays. 
Of these, the concentrations, in picograms per milliliter, of 
29 cytokines were measurable. The concentrations of the 
rest of the 19 cytokines were below the detection limit. 
Mean comparison procedures show that serum levels of 
six biomarkers (cutaneous T-cell attracting chemokine 
[CTACK], eotaxin, hepatocyte growth factor [HGF], 
interleukin 6 [IL-6], IL-16, and stem cell factor [SCF]; 
P<0.05) are significantly higher in COPD patients com- 
pared with in control subjects (Figure 1). The data of 
cytokines not showing significant differences between the 
groups and the rest of the biomarkers tested in the study 
are listed in Table SI. 
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Figure I Serum concentrations of cutaneous T-cell attracting chemokine, eotaxin, hepatocyte growth factor, interleukin 6 (IL-6), IL-I6, and stem cell factor are elevated in 
chronic obstructive pulmonary disease patients compared with control subjects. Serum level of (A) cutaneous T-cell attracting chemokine (P<0.0l), (B) eotaxin (P<0.05), 
(C) hepatocyte growth factor (P<0.05), (D) IL-6 (P<0.05), (E) IL-I6 (P<0.0l), and (F) stem cell factor (P<0.05) are significantly higher in chronic obstructive pulmonary 
disease patients (n— 2 1 ) compared with control subjects (n-9). Significant differences between chronic obstructive pulmonary disease and control subjects are denoted by 
*P<0.05 and **P<0.0l, as measured by two-tailed Mann-Whitney test. The data are represented as mean ± standard error. 
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CTACK is a T-cell chemokine, whereas eotaxin is a 
potent chemoattractant for eosinophils and basophils. 21,22 
HGF is a multifunctional mitogen. 23 IL-6 is a permeability- 
increasing inflammatory cytokine released mainly by 
macrophages and epithelial cells, and IL-16 is a proin- 
flammatory cytokine released by bronchial epithelium 
and dendritic cells, cluster of differentiation (CD)8 + and 
CD4 + T cells. 24 25 Hematopoietic factor SCF is secreted by 
various stroma and inflammatory cells. 26 Collectively, all 
these cytokines are important for recruiting different types 
of immune cells, cell proliferation, and tissue repair. The 
cytokines/biomarkers identified by these analyses are highly 
relevant to COPD and act as putative candidate biomarkers 
of COPD patients. 

Eotaxin is negatively correlated 
with FEV^ whereas CXCLI is 
positively correlated with HGF 

To determine the key biomarkers affecting FEV 1 (per- 
centage predicted), peak V0 2 (predicted), smoking, 
dyspnea score, and BMI in COPD patients, bivariate cor- 
relation analysis was performed. All recruited subjects 
had complete demographic data, and as expected by the 
study design, COPD patients have lower FEV, (percent- 
age predicted) (mean + standard error: 48.7%+3.9%) 
compared with age-matched controls (mean + standard 
error: 83.6%+1.7%; P<0.001; Table 1). Bivariate analysis 



Table I Demographic characteristics of the participants in the 
study* 



Parameters 


Healthy 


Patients with 


P-values 




controls 


COPD 




Sample size 


9 


21 




Age, years 


67.0 (S8.S-73.S) 


65.5 (53.0-76.0) 


0.09 


Forced expiratory 


83.6 (75.5-88.0) 


48.7 (32.5-67.0) 


0.001 


volume in 1 second, 








% predicted 








Peak oxygen uptake 




64.0 (44.0-81.0) 




(predicted) 








Dyspnea score 




2.0 (1.0-3.0) 




Body mass index 


22.0 (19.5-24.5) 


22.5 (19.5-23.5) 


0.15 


BODE-index 




5.0 (3.5-6.5) 




Smoking history or 


1 I.I (10.3-14.1)* 


18.0(14.0-23.5)' 


0.08 


status+ 









Notes: *AII the values are given as in median (interquartile range); *pack years; 
+only for smokers; ^healthy controls have no evidence of COPD (n-9): four and five 
subjects were smokers and nonsmokers, respectively. Patients with COPD (n-21): 
eleven of them are ex-smokers with COPD, whereas ten of them are current 
smokers with COPD. 



Abbreviations: COPD, chronic obstructive pulmonary disease; BODE, Body mass 
index, degree of airflow Obstruction and Dyspnea, and Exercise capacity. 
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Figure 2 Increasing severity of airflow limitation is associated with decreasing 
levels of eotaxin, whereas chemokine (C-X-C motif) ligand I (CXCLI) is positively 
correlated with hepatocyte growth factor (HGF). (A) Eotaxin is negatively 
correlated with forced expiratory volume in I second (FEy^ percentage predicted), 
as represented in the regression line. These data suggest that eotaxin is a good 
biomarker predicting FEV! changes in stable chronic obstructive pulmonary disease. 
(B) The regression line for CXCLI and hepatocyte growth factor shows a positive 
correlation, suggesting these two factors are dependent variables or regulated by 
the same transcription factor. 

of cytokines with each parameter shows that eotaxin 
concentrations have a strong inverse relationship with 
FEV, (percentage predicted; Y=0.106x + 87.08; r=-0.1l; 
P^O.05; Figure 2A). These data suggest that eotaxin could 
be one of the important contributors to the decaying lung 
function. Notably, chemokine (C-X-C motif) ligand 1 
(CXCL 1 ) or the growth-related oncogen a (GRO-a) show a 
strong positive correlation with HGF (Y=8.675x + 124.77; 
r=Q.16; P<0.05; Figure 2B); therefore, CXCLI and HGF 
are likely to be dependent variables and/or affected by 
the same transcription factor. CXCLI shows statistically 
significant but low levels of positive correlation with IL-6 
(r=0.23; P<0.05) and an inverse relationship with dyspnea 
(r=-0.53; P<0.05). CXCLI is a chemotactic factor for 
monocytes and neutrophils. This protein is produced by 
a variety of cells including monocytes, endothelial cells, 
and fibroblasts. 2728 Therefore, CXCLI may be involved in 
regulating various components of COPD. 
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Table 2 Multiple regression analysis between studied physio- 
logical parameters and cytokines showing significant differences 



between controls and patients with COPD 


Parameters and predictors 


Coefficient 


SE 


P-value* 


Forced expiratory volume in 1 second, % predicted 






Eotaxin 


-0.20 


0.02 


0.05 


Interleukin 6 


-0.37 


0.44 


0.04 


Stem cell factor 


-0.22 


0.03 


0.05 


Peak oxygen uptake (predicted) 








Eotaxin 


-0.46 


0.01 


0.04 


Cutaneous T-cell attracting 


-0.44 


0.02 


0.05 


chemokine 








Smoking history 








Eotaxin 


+0.45 


0.03 


0.04 


Hepatocyte growth factor 


-0.10 


0.01 


0.05 



Note: *A P-value ^0.05 is considered to represent a statistically significant factor. 
Abbreviations: COPD, chronic obstructive pulmonary disease; SE, standard error. 



Increase in CTACK, eotaxin, IL-6, 

and SCF significantly contributes 

to the reduction in FEN^ and peakV0 2 

To further determine the contribution of all the cytokines to 
each of these parameters, we conducted a stepwise multiple 
regression analysis. This systematic analysis identified that 
eotaxin, IL-6, and SCF are statistically significantly, contri- 
buting to the reduction in lung function measured by FEVj 
(percentage predicted). Contributions of each parameter are 
described by the following equation: Y=(-0.02xeotaxin) + 
(-0.37XIL-6) + (-0.22xSCF) (Table 2; P<0.05). A constant 
factor of — 1 .92 converts the Y values equivalent to observed 
FEVj (percentage predicted; Figure 3 A; r=-0.71; P<0.01). 
These analyses suggest that an increase in serum concentra- 
tions of CTACK, eotaxin, and SCF reflects the decrease in 
FEVj observed in the present study, and notably, these three 
cytokines are sufficient to explain most of the changes seen 
inFEV r 

Furthermore, eotaxin and CTACK significantly con- 
tribute to the observed reduction in peak V0 2 (predicted). 
Contributions of these two cytokines to peak V0 2 
are described by the following equation: Y=(-0.44x 
CTACK) + (-0.46xeotaxin) (Table 2; P<0.05). A constant 
factor of -9.1 converts the Y values equivalent to measured 
peak VO, (predicted; Figure 3B; ^-0.70; P<0.01). These 
analyses suggest that increase in CTACK and eotaxin con- 
centrations is related to decrease in peak V0 2 observed in 
these individuals. Once again, these two cytokines are suf- 
ficient to explain most of the changes observed in peak VO,. 
Collectively, these results show that the four key cytokines 
(CTACK, eotaxin, IL-6, and SCF) explain the significant 
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Figure 3 Cutaneous T-cell attracting chemokine, eotaxin, hepatocyte growth 
factor, interleukin 6 (IL-6), and stem cell factor concentrations estimate the 
changes in forced expiratory volume in I second (FEV ( ; percentage predicted), peak 
oxygen uptake (VO ; ; predicted), and smoking status by multiple regression analysis. 
(A) Stepwise multiple regression analysis identified that the variables eotaxin, IL-6, 
and stem cell factor significantly contribute to the changes in FEV (percentage 
predicted). The values calculated on the basis of these three cytokines correlate well 
with the changes observed in FEV! values. (B) Similar analysis shows that cutaneous 
T-cell attracting chemokine and eotaxin significantly contribute to the changes in the 
peak VO r The values calculated based on these two cytokines correlate well with 
the changes observed in peak V0 2 . (C) The smoking indices are related to eotaxin 
and hepatocyte growth factor, and the values correlate well with the predicted and 
observed smoking indices. 
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changes in lung function-related parameters observed in the 
study population. 

Smoking affects the changes in eotaxin 
and HGF concentrations 

To determine the effect of smoking on the changes in cytokine 
concentrations, we have conducted a stepwise multiple 
regression analysis of the biomarkers, as explained earlier. 
This analysis shows that changes in the sera concentrations of 
eotaxin and HGF are associated with the smoking levels (pack 
year): Y=(0.45xeotaxin) + (-0. IOxHGF) (Table 2; P<0.05). 
A constant factor of 12.6 converts the Y values equivalent to 
the observed smoking index (Figure 3C; /-0.58; 7 3 <0.05). 
These data suggest that smoking partially contributes to 
the increase in eotaxin and decrease in HGF concentrations 
observed in the sera of smokers in this study population. 

Specific immune and inflammatory 
response pathways are involved 
in stable COPD 

Complex interactions among cytokines, chemokines, 
and growth factors regulate various immune and inflam- 
matory responses. To determine the pathways that are 
potentially regulated by these cytokines in COPD, we 
have conducted enriched pathway analyses. Of the 
29 measurable biomarkers, 1 6 showed higher concentrations 
(£1.5 upregulated), whereas two showed lower concentra- 
tions (downregulated, sU.5) in COPD patients compared 
with in healthy controls. Pathway enrichment analysis was 
conducted using these cytokines by Metacore™ platform, 
with a priori setting of a cutoff threshold (False Discovery Rate 
(FDR)p, probability value =0.05). Pathways and biological pro- 
cesses regulated by these cytokines include immune responses 
(mediated by T-helper cell differentiation, riggering receptor 
expressed on myeloid cells 1 (TREM1) signaling, prostaglandin 
E2 signaling, andThl7 [T helper 17 cells] -derived cytokines), 
inflammation mediated by JAK/STAT (Janus kinase/signal 
transducers and activators of transcription) pathway, histamine 
signaling), lymphocyte proliferation and regulation, and leuko- 
cyte chemotaxis (Figure 4). Therefore, this analysis suggests the 
involvement of immune response-related pathways in COPD. 

GCR-a, NF-k(3 and C/EBP-p regulate 
the major cytokine networks in COPD 

Specific transcription factors and regulators determine 
cytokine response. 29 To determine the transcription 



factors that regulate the key cytokines involved in COPD, 
we have conducted a network analysis. Thirteen of the 
18 cytokines (nodes) generated a large single regula- 
tory network. Three main hubs NF-k(3 (nuclear factor 
kappa-light-chain-enhancer of activated B cells), C/EBP-P 
(CCAAT/enhancer-binding protein beta), and GCR-a 
(glucocorticoid receptor alpha) regulate these networks 
(Figure 5). These enriched networks mainly affect media- 
tors involved in immune, inflammatory, and chemotaxis 
responses. Interconnecting cytokine biomarkers including 
eotaxin, HGF, IL-6 and SCF are induced either by two 
or three hubs. Notably, GCR-a directly regulates all four of 
these cytokines, which are significantly upregulated in stable 
COPD patients. Therefore, GCR-a is likely a major potential 
regulator of these cytokine biomarkers in COPD. 

Discussion 

Serum biomarkers and their effects on COPD-related para- 
meters and pathways that affect COPD have not been fully 
established. We have assessed the levels of 48 serum bio- 
markers in COPD patients of Indian descent, using multiplex 
immunobead-based assays. We have identified some of the 
key biomarkers involved in pathways associated with stable 
COPD. In particular, the serum levels of CTACK, eotaxin, 
HGF, IL-6, IL- 1 6 and SCF are significantly higher in patients 
with COPD compared with in control subjects (Figure 1). 
These cytokine levels are expected to reflect long-term effects 
because our study focuses on stable COPD patients with no 
exacerbation during the previous 1 -month period. To the 
best of our knowledge, this is the first study that identified 
an increased serum level of SCF in COPD patients compared 
with healthy controls. Eotaxin, IL-6, and SCF inversely 
contribute to the changes in FEVj (percentage predicted), 
whereas eotaxin and HGF inversely contribute to the changes 
in peak V0 2 (predicted). Smoking is related to increased 
eotaxin and decreased HGF concentrations. CXCL1 is 
positively correlated with IL-6 and HGF, but it is negatively 
related to dyspnea. Enrichment pathways and network analy- 
ses revealed the involvement of inflammatory, chemotaxis, 
and immune response-related pathways in COPD. 

Previously, studies have reported that TNFa can contribute 
to COPD pathogenesis. Our data show a trend toward increased 
levels of TNFa in COPD patients (21.1-112.6 pg/mL) 
compared with in controls (21.1-34.1 pg/mL), but the 
values are not statistically different (Table SI). This could 
be a result of a lower level of TNFa in stable COPD patients 
or because of lower statistical power to detect the change in 
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Figure 4 Pathway enrichment analysis reveals specific immune and inflammatory response pathways in chronic obstructive pulmonary disease. (A) The top most significant enriched 
biological processes and pathways including immune responses, inflammation, lymphocyte proliferation and regulation, and leukocyte chemotaxis are represented along with the 
coenrichment P-value (log scale) in the bar graph. (B)The numbers and names of gene sets overlapping with the significant pathways and processes are shown in the table format. 
Abbreviations: PGE2, prostaglandin E2; IL-2, interleukin 2; IFN, interferon; GM-CSF, granulocyte-macrophage colony-stimulating factor; CCL2, chemokine ligand 2; 
SOS I, son of sevenless homologue I; VEGF, vascular endothelial growth factor; MIP-I, macrophage inflammatory protein I; HGF, hepatocyte growth factor; TREMI, 
triggering receptor expressed on myeloid cells I ; MIG, monokine induced by gamma-lnterferon; IP 1 0, interferon gamma-induced protein 1 0; MGF, mechano growth factor; 
STAT, signal transducers and activators of transcription; G-CSF, granulocyte colony-stimulating factor; Th 1 7, T helper I 7. 



this cytokine. Regulatory network analysis showed overlapping 
connectivity among eotaxin, HGF, IL-6, and SCF via 
C/EBP-(3, NF-k(3, and GCR-a transcription factors. In par- 
ticular, GCR-a is involved in the regulation of all four of 
these key cytokines in stable COPD patients. The involvement 
of these cytokines (eotaxin, HGF, IL-6 and SCF) in various 



pathways and networks indicates their potential regulatory 
roles in COPD. 

CTACK belongs to the C-C motif chemokines family, 
which plays an important role in the immune-inflammatory 
processes in many skin diseases. 21 CTACK is constitutively 
produced by epidermal keratinocytes and participates in 
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Figure 5 Interpathway interactions identify cytokine linking network in COPD. Diagram illustrates the network topology, in which individual annotated hubs (GCR-cc, 
nuclear factor kappa-light-chain-enhancer of activated B cells, and C/EBP-P) regulate the thirteen major cytokines. Note that eotaxin, hepatocyte growth factor, interleukin 
(IL)-6 and stem cell factor interconnecting nodes are induced either by two or three hubs. Red-dotted circles flag the upregulation, whereas blue dotted circles represent 
the downregulation of these cytokines in COPD. Green arrows show the positive effects, whereas red arrows represent negative effects. 

Abbreviations: COPD, chronic obstructive pulmonary disease; GCR-alpha, glucocorticoid receptor alpha; NF-kB, nuclear factor kappa-light-chain-enhancer of activated B 
cells; C/EBP-Beta, CCAAT/enhancer-binding protein beta chemokine (C-C motif) ligand 13; HGF receptor, hepatocyte growth factor receptor, met proto oncogene; MIP- la, 
macrophage inflammatory protein I alpha (CCL3); IL, interleukin; SOS I , son of sevenless homolog I ; MCP- 1 , monocyte chemoattractant protein- 1 ; G-CSF, colony stimulating 
factor (granulocyte); RANTES, regulated on activation normal T cell expressed and secreted; VEGF-A, vascular endothelial growth factor; MIG, monokine induced by gamma 
interferon (CXCL9); SCF, stem cell factor; GM-CSF, granulocyte-macrophage colony-stimulating factor; GOLD criteria, Global Initiative for Chronic Obstructive Lung 
Disease. 



tissue-specific homing of lymphocytes. 30 Increased level 
of CTACK present in COPD participants compared with 
controls (Figure 1A) suggests the involvement of CTACK 
in inflammatory processes related to the disease. In a study 
involving 48 patients, Pinto-Plata et al have also noted that 
CTACK levels are elevated in COPD subjects by protein 
microarray analysis. 31 Our study further shows that CTACK 
negatively affects peak V0 2 (Table 2). Therefore, including 
CTACK in future studies would be useful to confirm the 
relevance of this marker in other COPD populations. 

Eotaxin is a potent chemoattractant for human eosino- 
phils and basophils. 22 Our study shows that serum level of 
eotaxin in stable COPD patients is elevated (Figure IB). 
The COPD patients with eosinophil phenotypes show 
good response to corticosteroids. 32 D'Armiento et al have 



reported that eotaxin levels in plasma and bronchioalveolar 
lavage can predict rapid deterioration in lung function 
(FEVj) in COPD patients. 33 Similarly, our study shows 
an inverse relationship between serum eotaxin levels 
and FEV 1 (percentage predicted), as well as peak V0 2 
in the COPD subjects of Indian descent (Figure 3A and 
Table 2). Furthermore, smoking is positively corre- 
lated with eotaxin, suggesting that smoking potentially 
contributes to the increase of this cytokine. A previous 
study reports that smoking increases eotaxin levels in asth- 
matic subjects, 34 but the differences in eotaxin levels are not 
clearly established in COPD subjects. Pathway enrichment 
analysis further shows that transcription factors NF-k(3 and 
GCR-oc upregulate this cytokine production. Therefore, 
eotaxin is an important cytokine in COPD. 
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HGF is a multifunctional heterodimeric protein with 
mitogenic properties. 23 HGF is synthesized by fibroblasts, 
macrophages, smooth muscle cells, and epithelial cells, and it 
is an important cytokine for lung development and repair. 35 37 
HGF concentration is higher in COPD patients compared 
with in control subjects (Figure 1C). Smoking is negatively 
correlated with HGF (Figure 3C and Table 2), which could 
likely affect the repair process. Notably, the human HGF gene 
has an IL-6 response element, and therefore, HGF is expected 
to be upregulated during inflammation. 38 Collectively 
increased HGF identified in these patients may reflect the 
attempt by the host to repair the damaged lung. 

IL-6 is a classical marker for inflammation and is asso- 
ciated with many inflammatory disorders. 24 An increased 
concentration of IL-6 detected in stable COPD patients 
(Figure ID) is consistent with previous studies showing 
the increased level of this cytokine in the plasma of COPD 
patients. 39 Monocytes of COPD patients are known to over- 
react and release more IL-6 than cells from normal subjects. 40 
IL-6 is negatively related to FEVj (percentage predicted; 
Figure 3A and Table 2). Therefore, IL-6 is an important 
cytokine that could affect systemic inflammation and worsen 
COPD comorbidity. 4 

IL-16 is a proinflammatory biomarker released by bron- 
chial epithelium and dendritic cells, CD8 + and CD4 + T cells. It 
acts as a key element in COPD pathobiology. 25 IL-16 secreted 
by CD8 + cells is involved in regulating the recruitment and 
activity of CD4 + cells. The present study shows an increase 
in serum level of IL-16 in the COPD subjects (Figure IE), 
reflecting its potential involvement in the inflammatory 
process of the studied patients. Pinto-Plata et al have studied 
the serum biomarkers related to inflammation and injury in 
COPD patients and observed that the levels of the inflam- 
matory markers, including IL-6 and IL-16, were higher in 
more advanced disease. 41 The present finding is consistent 
with the involvement of these proinflammatory cytokines in 
the regulation of the disease process. 

SCF is a hematopoietic factor. 26 Various stroma and 
inflammatory cells express this factor. SCF stimulates mast 
cells and induces these cells to adhere to extracellular matrix, 
and lead to the production of proinflammatory cytokines 
and chemokines in the tissue. 42 - 43 For the first time, we have 
identified that COPD patients have higher levels of SCF than 
control subjects (Figure IF). This finding suggests that the 
SCF could directly remodel airway tissues in these patients. 
Dolgachev et al have shown that SCF, along with IL-31, is 
involved in remodeling and fibrosis of airways by promot- 
ing recruitment and differentiation of bone marrow-derived 



fibroblast precursors in mice sensitized and chronically 
challenged with allergens. 44 In patients with diffused inter- 
stitial lung fibrosis, alveolar fibroblasts secrete high levels of 
SCF. 45 COPD is also characterized by fibrosis and remodeling 
of airways; hence, SCF is likely to be a key factor involved 
in permanently altering lung architecture in COPD. 

CXCL1, or GRO-oc is a chemotactic factor produced by 
a variety of cells including monocytes, endothelial cells, and 
fibroblasts. Monocytes and neutrophils from patients with 
COPD show enhanced chemotactic response to CXCL1. 27 - 28 
Induced sputum of COPD patients has increased levels of 
CXCL1. 46 CXCL1 level shows negative correlation with 
HGF (Figure 3B) and positively correlates with IL-6 and 
dyspnea scores. The identified key cytokines/biomarkers 
may be involved in the regulation of the pathophysiology 
of the disease. 

The emerging field of network and pathways analy- 
sis can help identify the relationships between systemic 
inflammatory biomarkers and disease state. 47 - 48 The com- 
mon distinct functional and regulatory pathways identified 
in the present study are T H cell differentiation, TREM1 
signaling, prostaglandin E2 signaling, Thl7-derived cytok- 
ines, inflammation (mediated by the JAK-STAT pathway, 
histamine signaling), and leukocyte chemotaxis (Figure 4). 
One study shows that prostaglandin E2 is important for the 
mucosal innate immunity of COPD patients. 49 As evident 
from Figure 5, the network comprises three main seeding 
nodes (NF-k(3, C/EBP-(3 and GCR-a), which regulate 
inflammation and immune regulation. Many of the dif- 
ferentially expressed cytokines are regulated by any of 
these three nodes, and of them, IL-6, HGF, eotaxin and 
SCF are significantly higher in COPD subjects. A fuller 
understanding of these processes and pathways will pro- 
vide a better panel of COPD serum biomarkers and help 
us better understand the disease. 

One of the limitations of the study is the use of a relatively 
small sample size. Although 29 stable COPD patients who 
were selected with stringent exclusion criteria helped us 
identify key cytokines that are sufficient to explain most of 
the variability in FEV 1 and certain other parameters, and 
to identify pathways and networks, the cohort is not large 
enough to identify all of the potential biomarkers. This patient 
cohort is also not large enough to test comorbidities and to 
fully determine the effect of smoking on COPD biomarkers. 
Further studies with large sample size are necessary to 
validate these biomarkers, and molecular/cell biology-based 
studies are necessary to determine the mechanisms regulated 
by these cytokines in COPD. 



submit your manuscript | www.dovepress.com 
Dovepress 



International Journal of COPD 2014:9 



Dovepress 



Cytokine profiling and pathways analysis of COPD 



Conclusion 

In the present study, six biomarkers (CTACK, eotaxin, HGF, 
IL-6, IL-16 and SCF) are significantly elevated in the sera 
of COPD patients compared with healthy controls. Multiple 
regression analysis shows that eotaxin, IL-6, and SCF are 
inversely correlated with FEV r CTACK and eotaxin are 
negatively related to peak V0 2 . Smoking may increase eotaxin 
and decrease HGF Investigated enriched pathways and net- 
works suggest involvement of inflammatory JAK-STAT and 
histamine signaling, immune response medicated by TREM 1 
and prostaglandin E2 signaling and the regulation of lym- 
phocyte proliferation in the progression and pathogenesis of 
COPD. Notably, the four cytokines IL-6, HGF, eotaxin and 
SCF that are higher in COPD showed overlapping interactions 
with the seeding nodes. The involvement of these cytokines 
in various pathways and network analyses indicates their 
potential regulatory roles in the disease process. 
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Figure SI Flow diagram of the experimental design showing overall aim and analysis steps involved in the present study. 
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Table SI Concentrations (pg/mL) of biomarkers detected in the sera samples among groups 



Biomarkers (pg/mL) 


Non-smokers controls, n=S 


Smokers controls, n=4 


Ex-smokers with 
COPD, n=l 1 


Smokers with 
COPD, n=IO 


MIG 


907 


2,324 


5,201 


3,004 




(697-1,451) 


(1,874-3,174) 


(1,401-7,699) 


(1,738-5,413) 


IL-4 


3.4 


3.4 


4.0 


2.9 




(3.4-4.4) 


(3.1-4.0) 


(3.5-6.2) 


(1.3-4.1) 


IL-13 


10.1 


16.3 


16.6 


10.9 




(7.0-16.7) 


(13.8-21.8) 


(12.6-19.1) 


(8.2-15.0) 


SCGFfj 


17,831 


12,218 


26,283 


18,200 




(10,713-1,155) 


(8,122-21,971) 


(19,760-45,006) 


(12,617-3,655) 


IL-2 


7.3 


6.7 


1 1.8 


1 1.2 




(6.3-9.8) 


(4.8-8.5) 


(6.9-15.1) 


(8.4-14.5) 


GM-CSF 


24.1 


9.9 


49.8 


52.1 




(14.9-28.9) 


(4.3-14.6) 


(8.5-80.2) 


(16.9-92.5) 


RANTES 


4,226 


4,301 


2,951 


4,495 




(2,844-10,521) 


(3,181-5,145) 


(2,384.0-4,241.0) 


(3,806-8,808) 


CXCLI orGRO-a 


47.0 


48.7 


92.3 


1 16.2 




(20.2-1 1 1.2) 


(46.4-50.2) 


(58.5-1 14.5) 


(66.1-149.9) 


IL-12 


31.2 


63.6 


74.1 


37.9 




(19.0-60.7) 


(53.1-75.3) 


(30.5-92.4) 


(26.5-51.2) 


VEGF 


47.5 


80.4 


148.6 


89.7 




(36.8-101.8) 


(18.2-192.0) 


(92.9-194.4) 


(41.4-126.6) 


MlP-la 


137.2 


122.0 


103.8 


84.4 




(92.7-159.3) 


(36.4-264.0) 


(74.0-197.1) 


(63.5-98.1) 


IL-9 


30.2 


36.7 


42.3 


70.6 




(20.1-63.9) 


(27.6-60.5) 


(31.5-78.1) 


(32.9-657.8) 


FGF-p 


19.1 


16.0 


24.9 


17.6 




(15.9-25.3) 


(13.4-17.8) 


(14.4-40.3) 


(1 1. 1-25.1) 


MCP-I 


22.2 


24.6 


26.7 


33.3 




(15.3-24.3) 


(8.6-37.6) 


(15.5-41.9) 


(18.2^5.9) 


TNF-a 


25.7 


23.9 


32.7 


35.1 




(22.6-34.1) 


(21.1-29.9) 


(22.2-89.2) 


(21.6-1 12.7) 


PDGF 


10,103 


8,492 


10,862 


8,21 1 




(7,873-1 1,145) 


(7,774-10,431) 


(7,544-13,615) 


(5,795-10,298) 


IL-17 


82.2 


85.6 


96.7 


80.7 




(81.2-86.5) 


(79.2-95.6) 


(80.9-106.2) 


(67.7-93.4) 


IL-8 


9.4 


20.4 


14.9 


1 1.7 




(7.9^5.0) 


(7.3-191.4) 


(9.2-180.2) 


(9.9-20.3) 


IFN-Y 


178.1 


184.6 


206.6 


191.2 




(153.9-219.5) 


(15 I.I -208.5) 


(143.1-360.8) 


(109.7-246.5) 


IP- 10 


589.7 


722.4 


604.8 


565.4 




(353.5-734.1) 


(401.1-1,136.0) 


(268.1-1,488.0) 


(209.8-787.8) 


G-CSF 


14.1 


15.5 


13.7 


16.0 




(1 1.9-18.2) 


(1 1.5-17.3) 


(10.6-17.2) 


(1 1.3-19.6) 


TRAIL 


58.2 


42.7 


40.9 


70.8 




(9.3-88.4) 


(33.6-94.8) 


(27.8-86.2) 


(23.1-108.5) 


MIP-lp 


7.2 


7.5 


6.4 


6.6 




(5.2-9.1) 


(5.4-13.3) 


(5.8-14.5) 


(5.7-7.4) 



Notes: 'All the values are given as median (interquartile range). *These detected cytokines in the sera samples showed no significant differences in between either of the sub- 
groups including controls smoker and non-smoker and also in between ex-smoker with COPD and smoker with COPD. The following cytokines are undetectable {below the 
detection limit of the assay) in the sera: NGF-p, IFN-a2, IL-10, IL-15, IL-18, IL-la, IL-lp\ IL-IRa, IL-2Raa, IL-3, IL-5, LIF, MCP-3, M-CSF, MIF, MIP-lp\ SDF-la and TNF-p. 
Abbreviations: COPD, chronic obstructive pulmonary disease; MIG, monokine induced by gamma-interferon; SCGF, stem cell growth factor; GM-CSF, granulocyte 
macrophage colony-stimulating factor; RANTES, regulated on activation, normal T cell expressed and secreted; GRO, growth- related oncogene; VEGF, vascular endothelial 
growth factor; FGF, fibroblast growth factor; PDGF, platelet-derived growth factor; IP- 10, interferon gamma-induced protein 10; G-CSF, granulocyte colony-stimulating 
factor; NGF, nerve growth factor; IFN, interferon; IL, interleukin; R, receptor; Ra, receptor alpha; LIF, leukemia inhibitory factor; MCP, monocyte chemoattractant protein; 
M-CSF, macrophage colony-stimulating factor; MIF, migration inhibitory factor; MIP, macrophage inflammatory protein; SDF, stromal cell-derived factor; TNF, tumor necrosis 
factor. 
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